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We propose a high efficiency tomographic scheme to reconstruct an unknown quantum state of the qubits 
by using a series of quantum nondemolition (QND) measurements. The proposed QND measurements of the 
qubits are implemented by probing the the stationary transmissions of the dispersively-coupled resonator. It 
is shown that only one kind of QND measurements is sufficient to determine all the diagonal elements of the 
density matrix of the detected quantum state. The remaining non-diagonal elements of the density matrix can 
be determined by other spectral measurements by beforehand transferring them to the diagonal locations using 
a series of unitary operations. Compared with the pervious tomographic reconstructions based on the usual 
destructively projective (DP) measurements (wherein one kind of such measurements could only determine 
one diagonal element of the density matrix), the present approach exhibits significantly high efficiency for 
iV-qubit (N > 1). Specifically, our generic proposal is demonstrated by the experimental circuit-quantum- 
electrodynamics (circuit-QED) systems with a few Josephson charge qubits. 

PACS number(s): 42.50.Pq, 03.65.Wj, 03.67.Lx, 85.25.Cp 



I. INTRODUCTION 

The reconstruction of an unknown quantum state from 
a suitable set of measurements is called quantum tomogra- 
phy which is a particularly important method in the study 
of quantum mechanics and its various applications. Since the 
characterization of the states is the central tasks in quantum- 
state engineerings and controls, this technique is of great im- 
portance in the current quantum information processing. Re- 
cently, many theoretical analysis and experimental demonstra- 
tions have been devoted to implement the desirable quantum- 
state tomographies for, e.g., the polarization states of pho- 
tons miH, the electronic states of trapped ions 0], and the 
solid-state qubits |5], etc. However, all these tomographic re- 
constructions are based on the destructively-projective (DP) 
measurements, and thus are very operational-complicated. 
This is because that each kind of DP measurements, e.g., 
Pk = \k) (k\, is required to be performed many times on many 
copies of the reconstructed state for determining just one of 
the elements (e.g., |cfe| 2 in the state \ip) = J2k c k\k)) in the 
density matrix p = l^)^. 

Besides the usual DP measurements, quantum state could 
also be detected by other strategies, typically such as the so- 
called quantum nondemolition (QND) measurements. Basi- 
cally, QND measurement is a nondestructive detection, as the 
measurement-induced back-action noises could be effectively 
suppressed by repeatedly hiding them in certain observables 
which are not of interests. The basic criteria for a QND mea- 
surement is that the repeated measurements of an observable 
o of the same system should yield the identical result. This 
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means that the measured observable must be commutative 
with the Hamiltonian Hint describing the interaction between 
the measured system and the detector, i.e., [Hint, d] = 0. His- 
torically, QND measurement is proposed to explore the fun- 
damental limitations of measurements, and has been demon- 
strated in various branches of physics, such as in the detection 
of gravitational waves @], quantum optics telecommu- 
nications fioll . and quantum control ifllll . etc. In recent years, 
the QND measurement has also been successfully applied to 
probe the atomic qubits in the cavity quantum electrodynam- 
ics (QED) JT2L fl3fl . Furthermore, this technique was exten- 
sively used to the circuit QED systems IU4I420I1 for nonde- 
structively reading out the superconducting qubits. This QND 
measurement is implemented by measuring the transmission 
of the driven microwave signals through a transmission line 
resonator, which is dispersively coupled to the detected qubits. 
This is because that the detected qubits can cause sufficiently 
large state-dependent shifts of the resonator frequency. Thus 
by detecting the signals of the shifted frequency of the res- 
onator, the qubit state will be read out. However, the QND 
measurements in the above works lfl4l - l20ll are only utilized 
to effectively distinguish the different logic states of the de- 
tected qubit(s), which is (are) not prepared initially at their 
superposed state. 

Motivated by the above experiments, recently we proposed 
a new scheme to nondestructively detect the su perp osition of 
these logic states by the QND measurements 12111 . By tak- 
ing account of the full quantum correlations between the res- 
onator and dispersively-coupled qubit(s), our proposal shows 
that each detected peak marks one of the logic states and the 
relative height of such a peak is related to its corresponding 
superposed probability. This means that one kind of the QND 
measurements can determine all the diagonal elements of the 
density matrix of the measured quantum state p = 
Similarly, the non-diagonal elements of p could be determined 
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by other kinds of QND measurements by performing the suit- 
able unitary operations to transfer them into the measurable 
diagonal locations. Therefore, the proposed tomographic re- 
constructing approach is high efficient for N (N > 1) qubits, 
as the number of the kinds of the QND measurements required 
is significantly decreased. For example, to tomographically 
reconstruct a two-qubit state, the proposed 6-kind QND mea- 
surements are sufficient. This is significantly simpler than the 
previous schemes (requiring 15-kind measurements) based on 
either the D Ps pll or the individual dispersive readouts of the 
logic states ||20ll . 

The paper is organized as follows: Sec. II gives our generic 
model of the transmissions of a driven resonator. In Sec. Ill, 
we provide a detailed analysis of the QND measurement 
of a single-qubit state by probing the transmissions of the 
driven resonator. Next, we show how to use these QND mea- 
surements to tomographically reconstruct an unknown single- 
qubit state in the experimental circuit-QED system. The ex- 
tensions to the two-qubit case are given in Sec. IV, where the 
advantage of our proposal (compared with the previous ap- 
proach based on the DP measurements) will be explicitly re- 
vealed. The possible generalization to the 7V-qubit situation 
and summarizations of our main results are finally given in 
Sec. V. 



II. TRANSMISSION OF A DRIVEN EMPTY CAVITY 

For the detection of the states of the qubits, we investigate 
the photon transmission of a driven resonator by studying the 
steady-state properties of the resonator-qubits dynamics. For 
generality, we consider a cavity-QED system consisting of N 
qubits. The Hamiltonian reads 
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H = hu> r a f a + J2^^ L(T ^ + h 9i( a +i a + ^-j^)], (1) 

3=1 

where and a± j are ladder operators for the photon field 
and the jth qubit respectively. Also, tu r is the cavity fre- 
quency, ojj the jth qubit transition frequency, and gj the cou- 
pling strength between the jth qubit and the resonator. Sup- 
pose that the cavity is coherently driven by 
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where e is the real amplitude and uid the frequency of the ap- 
plied external driving. 

Under the Born-Markov approximation, the dynamics of 
the whole system with the dissipations and dephasings is de- 
scribed by the following master equation l22ll 
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Qn = —~[H N , g N ] + nV[a]g N + y^ j j ltj V[a_ j ]g N 

3=1 

N 

0) 

3=1 

Hn = H + Hd- 



Here, qn is the density operator and the dissipation super- 
operator is defined by 2? [A] g at = Ag^A^ — A^Aqn/2 — 
qnA'A/2, which describes the effects of the environment on 
the system. The parameters of the last three terms in Eq. (3) 
correspond to photon decay rate k, the jth qubit decay rate 
71 j, and the jth qubit pure dephasing rate 7^ j, respectively. 

In what follows, we begin with the master equation (3) to 
calculate the frequency-dependent transmission of the cavity, 
which is proportional to the steady-state mean photon number 
(era) in the cavity. Technically, to satisfy the basic criteria for 
the desirable QND measurements of the iV-qubit system, we 
assume that the conditions 
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«1, j^j' = l,2,...,N, 



(4) 

should be satisfied for assuring the effective dispersive cou- 
pling a Zj a^a between the jth qubit and the cavity. These con- 
ditions assure also that the inter-bit interactions are negligible. 
Above, Aj = ujj — u) r denotes the detuning between the jth 
qubit and the cavity, and Ajf = ujj—ujj the detuning between 
the jth and j'th qubits. 

For contrast, we first calculate the transmission spectrum 
of a driven empty cavity. The Hamiltonian of the simplified 
system reduces to (ft = 1 throughout the paper) 



H = Ur^a + eiflte**** + ae luldt ). 



(5) 



After the time-dependent unitary transformation defined by 
the operator R = cxp(— iutdta^a), we get the effective Hamil- 
tonian 

H = R^HqR - iR^dR/dt = -A dr tfa + e{a) + a), (6) 

where A f / r = ujd — w r is the detuning of the cavity from the 
driving. Consequently, we get the master equation for such a 
driven empty cavity 

g Q = -i[H , go] + n(ag a) - a)ag /2 - g a)a/2), (7) 

where go is the density matrix of the empty cavity. 

From the above master equation, one can easily obtain the 
equations of motion for the expectation values of the rele- 
vant operators, such as mean photon number inside the cavity 

(a^a) = Tr(a'aQo): 



d{a^a) 
dt 



K,(a'a) — 2elm(a), 



with 



d(a) ... 

The steady-state solution of Eq. (8) gives 

1 



(a t a) ss 



(8a) 



(8b) 



(9) 



Obviously, the transmission spectrum of an empty cavity, 
which is proportional to (a* a), is well-known Lorentzian: 
centered at u)d = oj r with the half-width k. Certainly, when 
u>d does not sufficiently match the cavity frequency, no photon 
penetrates the cavity and thus no transmission is recorded. 
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III. TOMOGRAPHIC RECONSTRUCTION OF A 
SINGLE-QUBIT STATE BY QND MEASUREMENTS 

A. Nondestructive detection of a single qubit by cavity 
transmissions 

Now we investigate the case, in which a single qubit with 
transition frequency uji is dispersively coupled to the cavity 
mode. In the frame rotating at drive frequency ui d character- 
ized by the transformation R, the Hamiltonian of the system 
reads 

Hi = yCT 2l + {-A dr + T 1 a Zl )a)d + e{a) +a), (10) 

with Ci>i = Wi + Ti and Ti = gf/Ai. 

Under the Born-Markov approximation, the master equa- 
tion for the single-qubit plus the driven resonator is 



Qi = -i[Hi, qi] + kD[o\qi + 7i. 1 £>[<T_ 1 ]£ 1 



(ID 



The desirable quantity (a^a) can be determined by solving the 
following coupled equations of motion: 



dt 



-Ac(a^a) — 2elm(a), 



(12a) 



Under the steady-state condition, we obtain 
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(A dr + T^a^ (0)))(KA dr + 7i,iA dr + 7i,i r i)] 



x[( 
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(nA dr 



A dr) 2 



7 M A dr + 7 M r 1 ) 2 ]-\ 
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which is strongly related to the the initial state of qubit, thus 
the qubit state could be determined by the cavity transmission. 




co d - co. (MHZ) 



d(d) k 

—— = (iA dr - -j){a) - iTi{da Zl ) - ie, (12b) 



d(da Zl ) 
dt 



and 



= (iA dr - - 



dt 



li,i){acr Zl ) - (iTi + 7i,i)(a) 
(12c) 



-Ti.iCKO + l)- 



(12d) 



It is obvious that the additional measurement-induced dephas- 
ing rate 7^1 does not influence the solution of the equa- 
tions. As the decay 71.1 of the qubit is significantly less 
than the decay rate K of the driven cavity, the average of a Zl 
could be safely assumed to be unchanged during the measure- 
ment. In fact, the characterized time of the detection is de- 
termined mainly by the decay of the cavity k. Such a quan- 
tity is about 2n x 1.69 MHZ |U9tj, which is obviously larger 
than 7i : i = 2ir x 0.02 MHZ [16]. Experimentally, the time 
interval of completing a single QND measurement is about 
T e = 40ns fl9ll , this is significantly shorter than the lifetime: 
T\ ~ 7.3//S and the decoherence time: T2 ~ 500ns lfl6ll . 
Therefore, during such a readout the decay of the qubit is neg- 
ligible, i.e., <<7 2l (T e )} = exp(-7i 1 iT e )(( t 7 zl (0)) + 1) - 1 w 
(^i(0)}. 



FIG. 1 : (Color online) Cavity transmission for the single-qubit states 
versus the probe detuning uid — u r . Five cases of the qubit states for 
|/3i | 2 = 0, 0.2, 0.4, 0.5, and 1 are shown. For comparison, the empty 
cavity (EMC) transmission is also plotted in black line. The peak 
shifts by — Ti or Fi correspond to single logic state |0) or |1). For the 
superposition states, the double-peak relative heights (in contrast to 
the peak height of the empty cavity transmission) present clearly the 
superposed probabilities of the two logic states. Here, the parameters 
are selected as: (ri,«,7i 1) = 2tt x (-7.38, 1.69, 0.02)MHz (H 

mi. 

The measured cavity transmission (normalized to the peak 
height of the empty cavity transmission) versus the probe 
frequency detuning are plotted in Fig. 1. Generally, the 
qubit is assumed to be prepared initially in the state \ip) = 
/3 |0) + /?i|l). Obviously, when /3 = 0, or 1, it reduces to 
the single logic state |1) or |0). Compared with the empty 
cavity transmission (plotted as the dark line in Fig. 1), one 
observes that qubit-resonator coupling leads to a right (left) 
shift of the single peak in the transmission spectrum by a 
quantity — T (T), which is dependent of the logic states for 
(cr Zl (0)) = -1 ((er 2l (0)) = 1). Thus, the shifts of the peaks 
can be used to mark the logic states of the qubit. However, 
when the qubit is in the superposition of the two logic states, 
e.g., |/?i| 2 = 0.2,0.4, and 0.5, respectively, we see that the 
situation is very different from the case for the single logic 
states. In this case, the spectrum shows two peaks whose lo- 
cations coincide with that for the single logic states, but the 
relative heights of these two peaks correspond clearly to the 
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superposed probabilities, i.e., \/3q\ 2 and |/?i| 2 , respectively. 
This provides an effective approach to directly measure the 
superposed probabilities of a superposed state. 



B. Tomographic reconstruction of a single-qubit state 

Above investigation indicates that, partial information of 
the qubit state, i.e., the diagonal elements of the relevant den- 
sity matrix, can be directly obtained by only one kind of the 
QND measurements. While, to extract the full information of 
an unknown qubit state, one should tomographically recon- 
struct all the elements of its density matrix. Basically, to com- 
pletely define a rf-dimensional density matrix p, one needs to 
determine d 2 — 1 real parameters. Therefore, to determine an 
unknown qubit state, the key point is to identify these param- 
eters by virtue of the tomographic technique. 

Now we demonstrate how to perform the tomographic con- 
struction of an arbitrary single-qubit state = /?o|0) + 
/3i 1 1), whose density matrix operator reads 

Pl = ( ^ J* ) . (14) 
\ Pw Pll J 

A more efficient and widely used technique is to parameterize 
the density matrix p\ on a Bloch sphere ||5|], 

i—x^y^z ' 

Here, / denotes the identity matrix, <jj the Pauli matrices, and 
r, real parameters. Therefore, in order to determine the single- 
qubit state, we must identify the three components (r x ,r y ,r z ) 
of the Bloch vector r. As discussed in the previous section, 
two diagonal elements poo an d Pn of the density matrix p\ 
can be directly determined by the two measured occupation 
probabilities |/?o| 2 , |/?i| 2 in the direct QND measurements. 
This means that the parameter r z can be determined by the re- 
lation r z = poo — Pn = /?o 1 2 — \Pi\ 2 - To obtain the other two 
parameters r x and r y , we need to determine the non-diagonal 
elements. To this end, we perform the single-qubit operations: 
U Xl = exp (ma Xl /A) and U yi = exp {iircr yi /&), to trans- 
fer them to the the relevant diagonal locations, respectively. 
For example, after the operation U Xl , the density matrix p\ is 
changed to 

Now, performing another kind of QND measurements the pa- 
rameters |/3 | 2 and \(3[\ 2 can be measured. Consequently, 
the coefficient r y can be determined via the relation r y = 
\f3'i\ 2 — \ftb\ 2 - Similarly, by performing the quantum operation 
U yi on the original density matrix py, another new density 
matrix 



can be obtained and the coefficient r x can be similarly deter- 
mined. Note that here the number of the unitary operations re- 
quired for implementing the quantum state tomography (based 
on the QND measurements) is the same as the previous ap- 
proach (based on the usual DP measurements). Thus, for the 
single-qubit case the complexity of the present approach is the 
same as that in the previous one. Note that here, as the same as 
that in the previous approach based on the DP measurements, 
three kinds of QND measurements are still required for the 
present reconstructions. One is directly applied, another is 
applied after the U Xl operation, and the final one is applied 
after the U yi operation, thus the efficiency is not enhanced. 

The remaining task is to implement the single-qubit oper- 
ations required above for transferring the non-diagonal ele- 
ments to the diagonal locations. We work with a circuit-QED 
system wherein a superconducting charge qubit is coupled to 
the fundamental mode of a transmission line resonator l23ll . 
Let the qubit work at its degeneracy point and neglect the 
fast oscillating terms under the rotating-wave approximation 
(RWA). Following Ref. l24ll . under one displacement transfor- 
mation, the effective Hamiltonian of the resonator plus qubit 
system can be written as 

H = —A dr a^a + -^-cr Zl + gi (a)a- 1 + aa +1 ) + — a Xl , 

(18) 

with the detuning of the qubit transition frequency from 
the drive A a = uji — w<j and the Rabi frequency ft = 
2egi/ (— Ad,-)- Next, supposing that this system works in the 
dispersive regime, i.e., |<7i/Ai| <c 1, after the transformation 
Ui = exp [-^(oV-j — aa +1 )/Ai], then the above Hamil- 
tonian becomes 

t A a tt 
H x = -A dr a 1 a+ + -^ a xi, A a = A a 

(19) 

First, if the condition A a = is satisfied, then the Hamil- 
tonian (19) produces a rotation of the qubit about the x axis, 
i.e., U Xl could be generated by choosing the evolution time 
t x = 7r/(2f2). Second, if the driving is sufficiently detuned 
from the qubit and its amplitude is also sufficiently large 
enough, then another approximate Hamiltonian 

1 1 f2 2 

H z = -A dr a)a+-{A a + - — )a Zl , (20) 

can be obtained by further performing a transformations 
U-2 = exp(/3*cr +1 — /3(T_ 1 ), with the coefficient /3 = 
n/(2A a ), on the Hamiltonian (19). Obviously, the de- 
sirable operation U Zl can be implemented by the evolu- 
tion under the Hamiltonian (20) with the duration t z = 
7rA a /(2A Q A a + f2 2 ). Third, the desirable operation 
U yi could be constructed as: U yi = cxp(i7rer yi /4) = 
exp(z7rcr Zl /4) exp(i37rcr a ; 1 /4) exp(i37rcr Zl /4). It should be 
pointed out that the durations t x (or t y ) of the single-qubit 
operations required above for implementing the desirable to- 
mographies is estimated as ~ lOOps using the experimental 
parameters: e - 2tt x 20MHz JH, and A dr ~ n/2 JH. This 
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is significantly less by at least two orders than the qubit deco- 
herence time, which is measured as ~ 500ns |24]. Therefore, 
the required gate operations are accessible and the proposed 
tomographic reconstructions are experimentally feasible. 

As an example, we assume that the three parameters r x = 
0.6, r y = 0.5, r z = 0.6 are obtained through the above recon- 
structions, then the reconstructed state p\ can be written as 
pi = 0.8|0}(0| + (0.3 - 0.25i)|0)(l| + (0.3 + 0.25i)|l)(0| + 
0.2 1 1) (1|, whose real p\j and imaginary p^' parts (i, j=0, 1) 
are graphically represented in Fig. 2. 




FIG. 2: (Color online) Graphic representations of the density matrix 



pi for a single-qubit state. The real p^' and imaginary p\'-' parts of 



the density matrix elements pij 
and (b), respectively. 



(i\p\j) (i, j=0, 1) are plotted in (a) 



and the equations of motion for the mean values of various 
expectable operators are 



d(a) 
dt 



d{acr Zl ) 
dt 



d(aa z 
dt 



d(a'a) 
dt 



-nla'a) — 2elm(a) 



(iA dr - -^)(d) - «Ti(aa Zl ) - iT 2 (da Z2 ) - it 



(23a) 
ie, 
(23b) 



(iA dr - 2 _ 7i,i)(oo- Zl ) - (iTi + ji,i){d) 



-iT 2 (aa zl a Z2 ) - ie(a zi ), 



(23c) 



(iA dr - - - 7i )2 ) (aa Z2 ) - (iT 2 +71,2) (6) 



-iT 1 (aa zl a Z2 ) - ie(a Z2 ), 



(23d) 



d(da Zl a Z2 ) 
dt 



(iA dr - - - 71,1 - 7i, 2 )(acr Zl cr Z2 ) 

-ie(<T Zl <T Z2 ) - (iT 2 + 7i,2)(oo"ai) 
-(*ri+7i,i)(a<r aa >, (23e) 



IV. TOMOGRAPHIC RECONSTRUCTION OF A 
TWO-QUBIT STATE BY QND MEASUREMENTS 



dt 



-n,i«f*i> + 1 ). 



(23f) 



A. Nondestructive detection of an unknown two-qubit state by 
cavity transmissions 

We extend the above sing-qubit QND measurements to the 
two-qubit case. The transition frequencies of the two qubits 
are represented as wi and u> 2 , respectively. In the above dis- 
persive condition (4) and in a framework rotating at uj d , the 
effective Hamiltonian of the present complete system is 

H 2 = (-A dr +T 1 a Zl +T 2 a Z2 )a' i a 



LO\ UJ'2 



w 2 



(21) 



where I\ 



j and uij 



Similarly, the relevant master equation reads 

Q2 = -i[H 2 , g 2 ] + KV[a]g 2 + 2J 7ij£ ) [c- j ]£»2 

J = l,2 



\Q2- 



(22) 



d(a z 



dt 



-71,2(^2) + 1), 



(23g) 



d(<J zl cr Z2 ) 
dt 



-(7i,i +7l,2)(o- 21 cr Z2 ) -7i,i(cr 22 ) 
-71,2 (a Zl ). (23h) 



Again, due to the relatively-long decoherence times of the 
qubits and their sufficiently short measured times, the addi- 
tional measurement-induced dephasing and decay rates of the 
qubits are also unimportant. Thus, the expectable values of 
the qubit operators can still be regarded as unchanged, i.e., 

K(*)> « K(°)> and K (*)"«*(*)} « (^(0)^(0)), 
during the desirable QND measurements. As a consequence, 
one can easily solve the above Eqs. (23a-e) and finally ob- 
tain the exact steady-state distribution of the intracavity pho- 
ton number 
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Here, A = iA dr - § , B 3 - = iT 3 , D 3 = iA dr Ej = 

iTj + 'fhjiF = i^dr - f - 71 J ~ 71,3'. G 3 = (<t Zj (0)), and 
Gi 2 = (a Zl (0)a Z2 (0))- 
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FIG. 3: (Color online) (a). Cavity transmission of the cavity 
versus the probe detuning uid — io r for certain sleeted two-qubit 
states with |a | 2 = 1, |ai| 2 = 1, \a 2 \ 2 = 1, a 3 | 2 = 1, 
(|a | 2 , |ai| 2 , |a 2 | 2 , |« 3 | 2 ) = (0.5,0,0,0.5), and (|a | 2 , |ai| 2 , 
|ct2 1 2 , |ct3 1 2 ) = (0.1,0.2,0.3,0.4), respectively. For compari- 
son, the empty cavity (EMC) transmission is also plotted in black 
line. Here, The parameters are {T\, Ty, k, 71,1, T1.2) = 2n x 
(-11.11, -9.11, 1.7,0.02,0.022, )MHz {JlH. In (b) only the 
parameters Ti and F2 are modified as r' x = 1.05Fi and T' 2 = 
0.85F2. In this case, the relative heights of the peaks are exactly 
equivalent to the corresponding probabilities of the single logic states 
superposed in the measured superposition state. 



We now investigate the above distributions schematically 
for various typically selected two-qubit initial states. First, 
we assume that the two-qubit is initially prepared at only 
one of the four logic states, i.e., in the generic expression 
\ip} 2 = ai|00)+a2|01)+a 3 |10)+a 4 |ll) only one of the four 
probability amplitudes equals 1, e.g., a\ = 1, a 2 = Oiz = &4 
= 0. Fig. 3 shows clearly that single peaks reveal the inputs 
of these four single logic states, and they can also be distin- 
guished by the shifts of the central frequencies of the transmis- 
sion spectrum. The peaks with frequency shifts: —Ti — T^, 
-Ti + T 2 , Ti - T 2 , and Ti + T 2 mark the state |00), |01), 
1 10) , and |11), respectively. Thus the pulls of the cavity are 
strongly dependent of the states of the qubits. For these sin- 
gle logic states the heights of the single peaks are exactly 
equivalent and of unity value, which is the same as that for 
the EMC. Next, for the superposition of the four single logic 
states the situations are quite different. For example, Fig. 3 (a) 



also shows that, if the two-qubit is prepared initially as one of 
the Bell states: (|ai| 2 , \a 2 \ 2 , \a 3 \ 2 , \a 4 \ 2 ) = (0.5,0,0,0.5), 
then the transmitted spectrum of the cavity reveals two peaks; 
their locations are respectively at the positions for the single 
states 1 00) and 1 11), but have the same relative heights. More- 
over, for a more generic superposed state (|ai| 2 , \a> 2 \ 2 , |cK3 1 2 , 
|a4| 2 ) = (0.1,0.2,0.3,0.4) one can see that four peaks are 
exhibited simultaneously. The central positions of these peaks 
locate at the corresponding positions of single logic states 
|00), 1 01 ) , 1 10) , and |11), respectively. The relative heights of 
them read 0.1, 0.212, 0.308, and 0.4, respectively. Here, the 
relative heights of the peaks marking the states |00) and 1 1 1 ) 
are exactly equivalent to the superposed probabilities |ai| 2 
and |ct4| 2 . However, the relative heights of the peaks marking 
the states 1 1 ) and 1 1 0) deviate from the corresponding super- 
posed probabilities \a 2 \ 2 and |ai3| 2 . This is because these two 
peaks are not well distinguished due to the contributions from 
these two logic states' overlap. As a consequence, each peak 
is higher a little than the expected one, i.e., the superposed 
probability of the relevant logic state. While, such a situation 
does not exist for the |00) and 1 1 1 ) peaks (the relative heights 
of them equal to the expected ones), as they are separated suf- 
ficiently far from the others. In Fig. 3 (b) we modify the rel- 
evant parameters such as T[ = 1.05Fi and T' 2 = 0.85F2. 
Then we find that each peak of the transmission of the cavity 
is well separated from the others, and thus its relative height is 
exactly equal to the expectable superposed probability of the 
corresponding logic state in the measured two-qubit state. 



High efficiency tomographic reconstruction of a two-qubit 
state 



The two-qubit state tomography is done in the same way 
as that for the single-qubit state. The only difference is that 
now there are 15 real parameters to be determined in the 4- 
dimensional density matrix operator p 2 , and thus more oper- 
ations are required to transfer the nondiagonal elements in p 2 
to the diagonal locations. Generally, the 4-dimensional den- 
sity matrix for a two-qubit state \ip) 2 = ai|00) + a2|01) + 
a 3 |10)+a 4 |ll) in a complete basis { 1 1)2 = |00), |2) 2 = |01), 
|3)2 = 1 10) , |4)2 = |H)} can be represented as 



Pi 



Pll Pl2 Pl3 Pli 

P21 P22 P23 P24 

P31 P32 P33 P34 

P41 P42 P43 P44 



(25) 



which can also be rewritten as p 2 ~ p 2 /4: with 
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p2 = 2_^i r mnO'm 1 Sv 0"„ 2 
m,n—0^x,y,z 

/ r 00 + r 0z + r z0 + r zz r 0x + r zx - ir 0y - ir zy r x0 + r xz - ir yQ - ir yz r xx - r yy - ir xy - ir yx \ 

rox + r zx + ir 0y + ir zy r 00 - r 0z + r z0 - r zz r xx + r yy + ir xy - ir yx r xQ - r xz - ir y0 + ir yz j 

r X Q + r xz + ir y0 + ir yz r xx + r yy - ir xy + ir yx r 00 + r 0z - r z0 - r zz r 0x - r zx - ir 0y + ir zy I 

V r xx - r yy + ir xy + ir yx r x0 - r xz + ir y0 - ir yz r 0x - r zx + ir 0y + ir zy r 00 - r 0z - r z0 + r zz ) 



Here, cr m =x,y,z are the Pauli operators and do identity matrix, 
and what we want to determine is sixteen real parameters r mn . 
Note that the first and second subscripts of the matrix elements 
Pij (i> j = l> 2, 3, 4) in Eq. (25) and r mn in Eq. (26) is labeled 
for the first and second qubits, respectively. 

As in the above discussion, performing the QND measure- 
ments on the two-qubit state p2 can directly determine all the 
four diagonal elements: pn, P22, P33 and ^44, respectively, 
by the measured results |«i| 2 , |«2| 2 , lo^l 2 an d | CK4 1 2 . As a 
consequence, the parameters vqq, r$ z , r z o and r zz can be de- 
termined by 



roo 
roz 

f ' ZZ 



\OL\\ 
l«l| 

|ai| 



F2I 
\a 2 \ 
\a 2 \ 
\a 2 \ 



\ota 



|a 4 | 
|a 4 | 

n 4 | 



1. 



using the similar approaches presented in Sec. III. With such 
a two-qubit operation and these single-qubit gates, we show 
in Table I how to perform the desirable unitary operations for 
transferring the non-diagonal elements to the diagonal loca- 
tions. For example, by performing a selected operational se- 
quence W = UffU Xi on the original density matrix p 2 , we 
have a new density matrix p' 2 = Wp 2 W', and the new diago- 
nal elements are 



(27) 



p'u = 


1 

4 (r °° 


1 Txy 


Tyz 


+ r zx ), 


p 22 = 


1, 

4 (r °° 


+ T x y 


+ r yz 




P33 = 


1, 

4 (r °° 


— Txy 


+ r yz 


+ r zx ), 


Paa 


\{rm 


— T xy 


— Tyz 





(30) 



To determine the other 12 parameters, we need to perform 
certain unitary operations to transfer them to the diagonal lo- 
cations for other QND measurements. 

It is well-known that, arbitrary two-qubit operation assisted 
by arbitrary rotations of the single qubits generate an univer- 
sal set of quantum gates. So the key to implement the above 
required operations for tomographies is to realize a two-qubit 
gate. Again, for the experimental circuit QED system with 
two superconducting charge qubits, such a gate could be im- 
plemented by using the so-called FLICFORQ protocol l24ll . 
In fact, if the cavity is driven by two external fields satisfying 
the sideband matching condition: uid 2 — u>di — ^1 + ^2> men 
an effective Hamiltonian 



Then, by the QND measurements the values of |o:^j 2 , la^l 2 , 
|«3| 2 , and \a' 4 \ 2 are given directly. As a consequence, the 
desirable parameters roo, r xy , r yz and r zx can be obtained by 
the relations: 



FF 



u> r a/a 



, 9l ff 2 (A , 1 + A^ 

16 a; a:. 



(28) 



can be induced in a quadruply rotating framework. Here, 

A j = u i +2Q 1o'/ A jdj> ~ w r with fijy = 2g j e j > / (u djl -u r ), 

and A jdr = w j ~ u}d i" = 1 > 2 ' - ? ^ i'- 0bviousl y' the 

evolution under the above Hamiltonian with the duration, e.g., 
around ~ lOOps, for the experimental parameters l24ll . can 
produce a two-qubit operation: 

U FF = cxp[i7r(cr yi (g) a y2 + a zi ® a Z2 )/4]. (29) 

On the other hand, the typical single-qubit gates, e.g., 
U Xj , U y ., and U Zj (j=l, 2) can be relatively easy to produce 



roo = \<xi\ ■ 

Txy = I | 2 

Tyz = -Kl 
I ' |2 

r zx = \a 1 \ 



J |2 



= 1, 



'3 



'4 



HI 



I4I 2 



(31) 



Similarly, other non-diagonal elements can also be deter- 
mined. Clearly, here only six kinds of QND measurements 
are sufficient to tomographically reconstruct a two-qubit state. 
This is obviously simpler than the previous tomographies 
based on the usual DP measurements, wherein 15 kinds of 
measurements are probably required I5I I20I1 . Thus, the present 
tomographies is essentially high efficient. 

After performing all the QND measurements listed in the 
table, a two-qubit state can be completely reconstructed. For 
example, a two-qubit state p 2 having the following represen- 
tation: 



Pi 



Too r 0x r 0y r 0z 

1 x0 r xx T xy T xz 

r v o r yx r yy r yz 

' z0 r zx T zy T zz 



( 1 



V 






-i 



K32) 



0/ 

which can be effectively reconstructed by these parameters 
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/ pn P12 pia Pu\ ( 0.1 0.0313 - 0.0313z 0.15 - 0.15* -0.125i 

Ptx P22 P23 P24 I _ [ 0.0313 + 0.0313* 0.2 0.125 -0.15 + 0.1 

psi P32 P33 P34 ~ 0.15 + 0.15* 0.125 0.3 -0.0313 + 0.0 

\P4i P42 P43 P44 / V 0.125* -0.15 - 0.15i -0.0313 + 0.0313i 0.4 



TABLE I: The operational combinations before the QND measure- 
ments to determine the parameters for tomographic ally reconstruct- 
ing a two-qubit state. The subscript "1(2)" of U is labeled for the 
operation of qubit 1(2). 



quantum operation W 



determined parameters 



no 

U FF U X1 
UffU v1 
U FF U Z1 
U v1 U z1 Uff 

Uy 2 U Z2 UFF 



Too, r zz , r 0z , r z0 

^"00. ^xy* Tyzi f zx 
TQQ, Tyx-, ?°zyi Txz 
^"00i Txx> Tyy, T zz 

n>o, rox, r y o, r yx 
n>o, r x o, ro y , r xy 



determined by six kinds of QND measurements. The simu- 
lated reconstructions are eranhicallv shown in Fie. 4. where 



P, 
st 
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(a) 



(b) 




FIG. 4: (Color online) Schematic representations of the density ma- 
trix p2 for a two-qubit state. The real and imaginary py parts 
(i.j = 1,2,3,4) of the density matrix elements in the complete 
bases are plotted in (a) and (b), respectively. 



V. DISCUSSIONS AND CONCLUSIONS 

Generally, the quantum state tomographic constructions 
demonstrated above can be extended for N (with TV > 2) 
qubits in a straightforward manner. This is because that the 
proposed QND measurements can be directly applied to deter- 
mine all the diagonal elements of the arbitrary A^-qubit state; 
the individual superposed logic states can be inferred from the 
relevant positions of the measured peaks, and the probabilities 
of the corresponding computational bases superposed in the 
measured state could be extracted from the relative heights of 
the peaks (when they separate sufficiently from the others). 
Moreover, all the required operations for the tomographic re- 
constructions can be implemented from the universal set of 
the logic gates demonstrated. 

In summary, we have proposed a scheme to perform the 
quantum state tomographies by QND measurements. Differ- 
ing from the usual tomographies based on the DP measure- 
ments, here the QND measurements are utilized. Since all 
the diagonal elements of the density matrix of an unknown 
quantum state can be simultaneously determined by a single 
kind of QND measurements, the efficiency of the present to- 
mographic reconstruction is definitely better for more qubits. 
Specifically, our proposal is demonstrated with the current cir- 
cuit QED setup with a few charge qubits, and could be gener- 
alized to other systems, in principle. 
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